Percolation of Color Sources and the Equation of State of QGP in central Au-Au 

collisions at ^/snn =200 GeV 
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The Color String Percolation Model (CSPM) is used to determine the equation of state (EOS) 
of the Quark-Gluon Plasma (QGP) produced in central Au-Au collisions at ^snn = 200 A GeV 
using STAR data at RHIC. When the initial density of interacting colored strings exceeds the 
2D percolation threshold a cluster is formed, which defines the onset of color deconfinement. 
These interactions also produce fluctuations in the string tension which transforms the Schwinger 
particle (gluon) production mechanism into a maximum entropy thermal distribution analogous 
to QCD Hawking-Unruh radiation. The single string tension is determined by identifying the 
known value of the universal hadron limiting temperature T c = 167.7 ± 2.6 MeV with the CSPM 
temperature at the critical percolation threshold parameter £ c =1.2. At mid-rapidity the initial 
Bjorken energy density and the initial temperature determine the number of degrees of freedom 
consistent with the formation of a ~ 2+1 flavor QGP. An analytic expression for the equation 
of state, the sound velocity Cf (£) is obtained in CSPM. The CSPM C,(£) and the bulk thermo- 
dynamic values energy density e/T 4 and entropy density s/T 3 are in excellent agreement in the 
phase transition region with recent lattice QCD simulations (LQCD) by the HotQCD Collaboration. 
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PACS numbers: 12.38.Mh; 25.75.Nq 



All high energy soft multi-hadron interactions exhibit 
thermal patterns of abundances characterized by the 
same temperature, independent of the center of mass en- 
ergy (Hi- The hadron limiting temperatures were mea- 
sured by statistical thermal analyses that fit the data 
with a minimum of parameters [l[ . In heavy ion collisions 
it may be plausible that multiple parton interactions pro- 
duce a thermalized system. In (p, p), (p, p) and (p, 
A) collisions multiple parton interactions are not likely 
to thermalizc the system. In (e + e~) annihilation where 
thermal behavior is observed, the multi-parton mecha- 
nism may also be an inappropriate explanation 043- 

The determination of the EOS of hot, strongly inter- 
acting matter is one of the main challenges of strong in- 
teraction physics (HotQCD Collaboration) [5j. Recent 
LQCD calculations for the bulk thermodynamic observ- 
ables, e.g. pressure, energy density, entropy density and 
for the sound velocity have been reported [5j]. In the 
present work a percolation model coupled with hydrody- 
namics has been utilized to calculate these quantities and 
compare them with the LQCD results. 

The CSPM describes the initial collision of two heavy 
ions in terms of a dense pattern of interacting colored 
strings which start to overlap to form clusters 0, Q . The 
interactions of strings reduces the hadron multiplicity fi 
and increases the average transverse momentum (p|,) of 
these hadrons, so that the total transverse momentum 
is conserved . The observables fi, and (pf) are directly 
related to the field strength in the string and thus to the 
generating color. For a cluster of n strings 
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is the average transverse momentum squared of a single 
string. The CSPM model calculation for hadron multi- 
plicities and momentum spectra was found to be in excel- 
lent agreement with experiment 8j. Within the frame- 
work of clustering of color sources, the elliptic flow, vg, 
and the dependence of the nuclear modification factor 
on the azimuthal angle show reasonable agreement with 
the RHIC data @. The critical density of percolation 
is related to the effective critical temperature and thus 
percolation may be the way to achieve deconfinement in 
the heavy ion collisions 10]. An additional important 
check of this interacting string approach was provided 
by the measurement of Long Range forward backward 
multiplicity Correlations (LRC) by the STAR group at 

RHIC (nj. 

The two dimensional (2D) percolation threshold iden- 
tifies the hadron to quark-gluon percolation phase tran- 
sition and the associated critical temperature. The 
Schwinger barrier penetration mechanism modified by 
Gaussian fluctuations about the average value of the 
string tension leads to an early (at birth) thermal system 
[Tol . Above the percolation threshold this system is 
considered to be in the deconfined phase, which subse- 
quently expands according to Bjorken boost invariant ID 
hydrodynamics |13| . 

With an increasing number of strings n there is a pro- 
gression from isolated individual strings to clusters and 
then to a large cluster which suddenly spans the area. In 
two dimensional percolation theory the relevant quantity 
is the dimensionless percolation density parameter given 

by@,a 



where //q is the multiplicity of a single string and (p£) 
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where N is the number of strings formed in the collisions 
and Si is the transverse area of the a single string and 
Sn is the transverse nuclear overlap area. The critical 
cluster which spans Sn, appears for £ c > 1.2 [TJj]. As 
£ increases the fraction of Sn covered by this spanning 
cluster increases. 

The color suppression factor F(£), reduces the hadron 
multiplicity from n/zo to the interacting string value [i 
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To evaluate the initial value of £ from data, a param- 
eterization of p-p events at 200 GeV is used to compute 
the p t distribution fl5| 



dN c /dp 2 t =a/(p +PtY 
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where a, poi and Qf are parameters used to fit the data. 
This parameterization also can be used for nucleus- 
nucleus collisions to take into account the interactions 
of the strings fl5| 
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where S n corresponds to the area occupied by the n over- 
lapping strings. The thermodynamic limit, i.e. letting n 
and S n — > co while keeping £ fixed, is used to evaluate 
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In pp collisions (nS\/ S n ) pp ~ 1 due to the low string 
overlap probability. The factor 1 — in Eq. 4 cor- 
responds to the fractional area covered by the spanning 
cluster. The STAR analysis of charged hadrons for 0- 
10% central Au+Au collisions at ^/snn =200 GeV gives 
a value £ = 2.88 ± 0.09 [JJ]. 

The strong longitudinal chromo-electric fields produce 
Schwinger-Bialas [l?J, [H, [3 like radiation with a ther- 
mal spectrum, in analogy with the Hawking-Unruh radi- 
ation [l7l - [22l |. Both the Schwinger-Bialas and Hawking- 
Unruh derivations lead to the same value of the maxi- 
mum entropy temperature. Above the critical value of 
£, the QGP in CSPM consists of massless constituents 
(gluons). The percolation parameter £ determines the 
cluster size distribution, the temperature T(£) and the 
transverse momentum in the collision The connec- 
tion between £ and the temperature T(£) involves the 
Schwinger mechanism (SM) for particle production. In 



CSPM the Schwinger distribution for massless particles 
is expressed in terms of p 2 
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with the average string value of (x 2 ). Gaussian fluctua- 
tions in the string tension ( Bialas) around its mean value 
transforms SM into the thermal distribution 1121 
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with (x 2 ) = ir (pi) i / 'F '(£) . The temperature is given by 
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In the determination of temperature using Eq.(ll) the 
value of F(£) is obtained using the experimental data 
[l5| . We will adopt the point of view that the experi- 
mentally determined chemical freeze-out temperature is 
a good measure of the phase transition temperature, T c 
[23j. The single string average transverse momentum 
(p 2 )i is calculated at £ c = 1.2 with the universal chemi- 
cal freeze-out temperature of 167.7 ± 2.6 MeV [l|. This 
gives v /(p2) 1 = 207.2 ± 3.3 MeV which is close to ~200 
MeV used previously in the calculation of percolation 
transition temperature [l0|. Above £ c =1.2 the size and 
density of the spanning cluster increases. We use the 
measured value of £ = 2.88 to determine the tempera- 
ture, before the expansion, 7$ = 193.6±3.0 MeV of the 
quark gluon plasma in reasonable agreement with Ti = 
221±19 stat ± 19 sys from the enhanced direct photon ex- 
periment measured by the PHENIX Collaboration [24j. 

The QGP according to CSPM is born in local thermal 
equilibrium because the temperature is determined at the 
string level. We use CSPM coupled to hydrodynamics to 
calculate energy density, pressure, entropy density and 
sound velocity. As mentioned earlier the strings interact 
strongly to form clusters and produce the pressure at the 
early stages of the collisions, which is evident from the 
presence of elliptical flow in CSPM Q. After the initial 
temperature T > T c the CSPM perfect fluid may expand 
according to Bjorken boost invariant ID hydrodynamics 
131 
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s = (l + C 2 )^ 
de s 
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where e is the energy density, s the entropy density, t the 
proper time, and C s the sound velocity. 
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Above the critical temperature only massless particles 
are present in CSPM. The initial energy density £j above 
T c is given by [H] 



2 S n Tp ro 



(17) 



To evaluate £j we use the charged pion multiplicity 
dN c /dy at midrapidity and S n values from STAR for 
0-10% central Au-Au collisions with ^/snn =200 GeV 
[25j . The factor 3/2 in Eq.(17) accounts for the neu- 
tral pions. We can calculate (pt) using the CSPM ther- 
mal distribution Eqs.(10) and (11). For 0.2 < p t < 1.5, 
(p t ) = 0.394 ± 0.003GeV, adding the extra energy re- 
quired for the rest mass of pions at hadronization (nit) — 
0.42 ± 0.003GeV. The error on (p t ) is due to the error on 
Ti. 

The dynamics of massless particle production has been 
studied in QE2 quantum electrodynamics. QE2 can be 
scaled from electrodynamics to quantum chromodynam- 
ics using the ratio of the coupling constan ts I2q| . The 
production time T pro for a boson (gluon) is [27( 
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Using Eqs. (17) and (18) gives = 2.27±0.16 GeV//m 3 
at £=2.88. In CSPM the total transverse energy is pro- 
portional to £. From the measured value of £ and e 
it is found that e is proportional to £ for the range 
1.2 < £ < 2.88, ei = 0.788 £ GeV//m 3 [H H|. This 
relationship has been extrapolated to below £ = 1.2 and 
above £ = 2.88 for the energy and entropy density calcu- 
lations. The Au-Au at ^/sWn =200 GeV data is used to 
normalize the CSPM e/T 4 values. Figure 1 shows e/T 4 
as obtained from CSPM along with the LQCD calcula- 
tions [5] and the CSPM pressure 3p/T 4 . 

The number of degrees of freedom (DOF) are related 
to the energy density 
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At T t the DOF is 37.5±3.6. At T c , e c =0.95±0.07 
GeV/ /to 3 and 27.7±2.6 DOF. 

The sound velocity requires the evaluation of s and 
dT/de, which can be expressed in terms of £ and F(£). 
With q 1 / 2 = F(£) one obtains: 



dT _ dT dq d£ 
de dq d£ de 



Then C 2 becomes: 



C 2 = (l + C s 2 )(-0.25) 
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for £ > £ c , an analytic function of £ for the equation of 
state of the QGP for T > T c . 
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FIG. 1: The energy density from CSPM versus t/T? spm 
(red circles) and Lattice QCD energy density vs T/T^ QCD 
(blue dash line) for 2+1 flavor and p4 action [IJ. 3p/T 4 is 
also shown for CSPM with green circles. 
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FIG. 2: The speed of sound from CSPM versus T/T C CSPM (red 
circles) and Lattice QCD-p4 speed of sound versus 
T/T^ QCD (blue dash line)!!]. The physical hadron gas with 
resonance mass cut off M < 2.5 GeV is shown as solid green 
circles HJ. 



Figure 2 shows the comparison of C 2 from CSPM and 
LQCD. The LQCD values were obtained using the EOS 
of 2+1 flavor QCD at finite temperature with physi- 
cal strange quark mass and almost physical light quark 
masses Q. At T/T c =l the CSPM and LQCD agree with 
the C 2 value of the physical hadron gas with resonance 
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mass truncation M < 2.5 GeV [28j . 

The entropy density s/T 3 is computed using Eq. (15) 
as shown in Fig. 3 along with the LQCD results. CSPM is 
in excellent agreement with the LQCD calculations in the 
phase transition region for T/T c <1.5. It is noteworthy 
that CSPM at £ > £ c exhibits gluon saturation effects 
similar to the Color Glass Condensate (CGC) [29|. The 
saturation scale Q s in CGC corresponds to (pt)i/F(0 m 
CSPM. 
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FIG. 3: The entropy density from CSPM versus 
T/T° SPM (red circle) and Lattice QCD entropy density vs 
T/Tt QCD (blue dash line) [|. 

In summary, the present analysis consists of using the £ 
values from data, identifying the universal chemical tem- 



perature with the critical percolation density parameter 
£ c , establishing the direct proportionality between £j and 
£ and using the functional dependence of T(£) and e, to 
obtain C 2 S . Thus we can determine the temperature de- 
pendence of energy density, pressure and entropy density 
in the color string percolation picture. The results are 
also in agreement with lattice QCD in the phase transi- 
tion region, when the results are plotted with respect to 
T/T° SPM and T/T^ CD . The value of C s 2 =0.14 is in 
agreement with the physical resonance gas value at the 
critical temperature. The non-interacting high tempera- 
ture limit for C 2 S = 0.33 is reached at T ~ 1.5T C . 

The percolation critical transition is known to repre- 
sent a continuous phase transition. In central Au-Au col- 
lisions at y/spjN = 200 GeV the QCD to hadron phase 
transition for baryon density /is ^0 is believed to be a 
cross-over transition which does not have a latent heat 
[H. The CSPM EOS correctly describes the QCD to 
hadron cross-over transition and provides an answer to 
the question of the origin of universal temperature ob- 
served in A-A, p-p and e + e~ collisions. 

The percolation analysis of the color sources applied 
to STAR data at RIFIC provides a compelling argument 
that the QGP is formed in central Au-Au collisions at 
y/sJ?N — 200 GeV. It also suggests that the QGP is pro- 
duced in all soft high energy high multiplicity collisions 
when the string density exceeds the percolation transi- 
tion. A further definitive test of CSPM can be made at 
LHC energies by comparing hadron-hadron and nucleus- 
nucleus collisions. 
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